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Abstract 

We study normalization problems associated with use of perturbatively correlated ground- 
states in extended RPA schemes in the context of a specific but typical example. The 
sensitivity of the results to the amount of 2p2h admixtures to the correlated ground state is 
also investigated in terms of a modification of the standard perturbative approach. 



PACS numbers: 21.60.Jz, 21.10.Re, 25.40.Ep 



Some time ago Van Neck et al.jJT] pointed out that some annoying numerical inconsis- 
tencies result from the evaluation of consistently derived perturbative expressions in the 
context of the nuclear many-body problem. Specifically, they pointed out that an often used 
procedure of evaluating the number of nucleons perturbatively excited above the Fermi level 
(including normalization factors expanded to the appropriate perturbative order) leads to 
grossly overestimated results. This happens due to the fact that, as the perturbation adds 
a very large number of relatively small excited 2p2h components to the OpOh wave function, 
the relative weight of the former in the perturbed wavefunction is typically large enough 
numerically so that the perturbatively expanded normalization becomes inadequate. This 
type of difficulty affects also linear response calculations done in the context of the so called 
extended second random phase approximation (ESRPA) 0, which uses a perturbatively gen- 
erated ground state wavefunction with 2p2h admixtures, in addition to including two-body 
operators in the structure of the phonons. In this note we work out an example that shows 
that this numerical normalization error tends in fact to inflate significantly ESRPA strength 
distributions, as was also pointed out in Ref. |IJ]. Moreover, considerable excess strength 
still remains over the results obtained by using just the unperturbed ground state (second 
random phase approximation (SRPA)) when one attends to the numerical normalization 
problem. This excess strength appears to be related to the relative importance of the 2p2h 
admixtures to the unperturbed ground state. 

In order to explore this point we give also results obtained for a modified ESRPA in which 
2p2h ground state correlations are introduced by means of the Brillouin-Wigner (BW) pertur- 
bation theory, which has the effect of reducing appreciably their importance. This happens 
through the lowering of the ground state energy produced by solving the appropriate disper- 
sion equation. The resulting strength distribution, for reduced but still non negligible 2p2h 
admixtures to the unperturbed ground state, comes out close to the simple SRPA result. We 
take these facts as an indication of enough sensitivity of the calculated strengths to the corre- 
lation structure of the ground-state so as to warrant the development and implementation in 
realistic situations of better controlled extensions of the standard quasi-boson random-phase 
approximation. 

We base our argument on the linear response R(E) to an external field F, which admits 
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the spectral representation 

R(E) = E 



(0|F|i/)(i/|Ft|0) _ (0\F^\u)(u\F\0) 
E — E y + irj E + Ey-if] 



Here |0) and \v) are the exact ground state and excited eigenstates of the full hamiltonian 
H. The excited states \v) can be written as 



(2) 



where the set {Ci,Cj} constitutes a complete operator basis. In the SRPA this set is re- 
stricted to one and two particle-hole annihilation and creation operators out of the Hartee- 
Fock (HF) ground state \HF) and the coefficients X" and Y? are determined from the 
equations of motion M 



(HF\ [ty,, \hM]\ \HF) = E U (HF\ fn M ,ntl \HF) = EJ^. 



(3) 



The ESRPA hinges on the idea that the inclusion of 2p2h operators among the Ci requires 
a modification of the quasi-boson approximation in which the HF ground state is allowed to 
include perturbative 2p2h correlations. Hence, in evaluating Eq. (|3|) one uses in this case a 
ground state of the form 



|6> = c 



|#F) + Ec2 |2o> 

2o 



(4) 



where the amplitudes C2 are evaluated in the first order Rayleigh-Schrodinger (RS) pertur- 
bation theory, i.e., 



C2 



(2 \V\HF) 



-E 



2 



(5) 



Here 2 = (^1^2^1^2)0 indicates 2p2h excitations with independent-particle energy E 2o , V is 
the residual interaction, and Cq is a normalization factor. Since Eq. (|4]) is thought of as a 
perturbatively generated expression, cq is generally set equal to 1. 
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In general Eq. (|3|) leads to a secular problem of the form 



AX" = E V NX 



V 



with 



A 



A B \ 

b* a* r 



X u 



x u \ 

yv , 



N 



N 
—N* 



(6) 



(7) 



where the submatrices A, B and N given by 



Aij — (0| Cj, 



a 



h,cX\ |o), a^ = (6| [a,ctl |0> 



H,C}\\\0), B id = (0\ 
Furthermore, one can write Eq. (H) in a representation independent form as 

R(E) = J* (EN + ir}2 - A)- l T, 
where the matrix T represents the operator F and is defined as 



(9) 



F A 

F B 



with 



if =<0| [Ci,F\ |0>, 

F B = F A*(p _^ p\y 



(10) 



Eq. (Bp can be reduced with the help of projection operators P and Q onto subspaces 
involving lplh and 2p2h excitations respectively. One gets 



R(E) = Fl{E)Q P {E)Fp{E) + F ] Q Q Q {E)F Q , 



where 



Q P {E) = [EN P + irfip -A P - (Apq - N PQ E) Qq(E) (A qp - NqpE)]' 1 , (12) 



with 



and 



G Q {E) = [EN Q + ir 1 l Q -A Q }-\ 



F P (E) =T P - Np Q T q + ApqQ q {E)T q . 



(13) 
(14) 



When using the ESRPA some more complicated two-body effects are trimmed by keeping 
terms up to second order in V for the forward sector within the P space, terms linear in V 
for the backward sector within the P space and for the coupling between the P and Q spaces, 
and only terms of zeroth order within the Q space. The usual argument (see e.g. Ref. H) 
for this procedure involves again the limitations stemming from the perturbative dressing of 
the ground state, Eq. (|). 

Finally, we set up a modified extended second RPA (MESRPA) in which Eq. (|J) is 
replaced by the corresponding expression obtained from the BW perturbation theory. This 
in fact coincides with the ESRPA result but the coefficients C2 are now given as 

W\HF) 

c 2o = -Ei p— , (15) 

^0 - ^2 

where the ground state energy E is the lowest solution of the secular equation 

E ^\Mwi)l, (16) 

2o ^0 - -^2 



This leads to increased energy denominators in Eq. (|15|) and hence to reduced 2p2h admix- 
tures to the HF ground state. One obtains in this way 

Nij = Sij + ANij (17) 

where % = ipih and the nonzero ANy are just 

AN 1V = \c \ 2 ^<2 |I>ii'|2o>, (18) 

where D\y = Ci,C\, — 5iy. (Note that within the quasi-boson approximation Diy = 0). 
The explicit result for the matrix element (2q|-Dii'|2q) is 

((piPhhUD^Mp^h^o) = - [1 + P(h x , h 2 )P{h v M')\ (19) 
x Sp,p'hiwP~(h,h2)P~(pi,P2)Sh 1 ,,h5h 2 ,h 2 ,8 P2 ,p 2l S PuPl , +p^h, 
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where P = [1 — P(i,j)], while the operator P(i,j) exchanges the arguments i and j. 

The matrix elements of A are 

A i:j — dijEj + Vy + AAij, (20) 
where V{j = (iV\j) and the nonzero matrix elements AA^ are: 

AA 1V = |c | 2 (Ei - E 2 + E )c; o% (2 Q \b n ,\2' ). (21) 

2o,2' 

Finally the matrix elements of T are: 

| A + Ei' ANii'fv for 2 = 1 (22) 
\ CoE 2o c 2o/22o for i = 2, 

(l\F\HF), and / 22o = (2|F|2 ). (23) 

Note that the corresponding ESRPA quantities are obtained by setting cq = 1 and -E^ = 
in the MESRPA expressions. 

We next give numerical results for the GT resonance (F± = at±) in 48 C a using the MY3 
force in the Oku — 3hu oscillator space. Four different ways of handling the nuclear ground 
state will be compared. The first one is just the plain SRPA in which the equations of motion 
(H) are evaluated with the HF ground state. We give also results for the ESRPA (for which the 
normalization coefficient Co is set equal to 1), for a normalized version of this approximation 
(NESRPA) in which cq is determined so that (0|0) = 1 with the C2 coefficients given by Eq. 
(P), and finally for the MESRPA, which uses a normalized ground state with amplitudes C2 
evaluated using Eq. dl5|). In order to obtain smooth strength functions S(E) = —-Im R(E) 
with R(E) given by Eq. (|Tl|), the energy variable is taken to be complex: E —>■ E + iA, 
with A = 1 MeV for the lplh and A = 3 MeV for the 2p2h subspace respectively. Solving 
the dispersion equation Eq. ([16|) gives E = —29 Mev, which amounts to about 8% of the 
experimental ground state binding energy. The results are shown in Fig. 1 and in Table 1 
below. The positive branch of the GT sum rule S + — <S_ = 3(iV — Z) with S± = J S±(E)dE 
is divided into a low energy part 5^ (E < 20 Mev) and a high energy part (E > 20 Mev). 



Ft 

where 
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These quantities are the relevant ones for the problem of the quenching of GT strength. In 
the usual RPA the low energy part 5+ essentially exhausts the sum rule. When lplh — 2p2h 
coupling is introduced via the SRPA, 30% of this strength is shifted to the high energy 
region. This amount is somewhat reduced when ground state correlations are introduced 
via the ESRPA. This results from the combined effect of the Q-space part of Eq. (|TT|) and 
of the interference effects generated by the last term of Eq. (|14]) || . Furthermore one gets 
now also a contribution in the negative branch S- so that increases to 82% (third line 
in Table 1). As shown in the first two columns of Table 1, the ground state wavefunction 
involved in the derivation of the ESRPA expressions has a serious normalization problem. 
In the NESRPA this is fixed by suitably reducing the value of cq. This has only a relatively 
small effect on <S+ and reduces both and the negative branch contribution 5_ (fourth 
line of Table 1). When ground state correlations are introduced via the MESRPA, on the 
other hand, the percentage of 2p2h admixtures to the ground state is reduced from 62% to 
31% while the strength distribution becomes quite similar to the simple SRPA result. This 
last feature indicates important sensitivity to the amount of ground state correlations which 
therefore, as stated above, deserves a more controlled treatment. It is worth noticing that, 
when the interaction among the 2p2h configurations is neglected, the BW approximation for 
the ground state wavefunction coincides with the diagonalization procedure |5|. This means 
that the OpOh — 2p2h coupling in the initial nucleus is treated at the same footing as the 
lplh — 2p2h coupling in the final nucleus. Thus we feel that, in the context of the present 
calculations, it is more consistent to use the BW perturbation theory than the RS one. 

Even though the present discussion has been limited to one specific case involving the 
Gamow- Teller response within the extended RPA, the observed trends should apply also to 
other schemes which include ground state correlation effects perturbatively, both for this |7j] 
and for other types of response functions, notably the longitudinal and transverse inclusive 
responses in quasi- free electron scattering || [H|]. In all cases ground-state normalization 
is numerically important and sensitivity to the amount of 2p2h correlations should be high, 
so that a moderate reduction of the 2p2h ground state component will lead to results which 
are not far from those obtained in the simple SRPA. 
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TABLES 



Table 1: Ground state normalization factor, summed weights of 2p2h components (see Eq. 
(£|)), GT integrated strength S + in the resonance region (<), above it (>), and total GT 
strengths S + and 5_. The first column identifies the approximation scheme. Strengths are 
given in percent of 3(N — Z). 
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RPA 


1 





100 





100 





SRPA 


1 





70 


30 


100 





ESRPA 
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1.60 


82 


27 


109 
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NESRPA 


0.38 


1.60 


80 


23.5 


103.5 


3.5 


MESRPA 


0.69 


0.45 


72.8 


28 


100.8 


0.8 



FIGURES 



Figure 1: Folded Gamow- Teller strength distributions in Ca for different approximation 
schemes. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9403017vl 



